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VESSELS PARTIALLY SUPPORTED BY SOIL 


W. A. Boothe,’ R. T. Gray,’ and G. Horvay’ 


ABSTRACT 


Curves are plotted which indicate the reduction in bending stresses that can 
be achieved near a discontinuity in a cylindrical shell by application of partial 
elastic support in the vicinity of the discontinuity. 


INTRODUCTION 


Cylindrical pressure containers occasionally are supported by encasing 
their bottom portion in concrete. In order to reduce the bending stresses that 
develop at the interface where built-in and free portions join, some manufac- 
turers have adopted the practice of introducing a transition zone of partial re- 
straint. Sand of suitable stiffness dammed up to a depth of several shell at- 
tenuation lengths forms a buffer layer, and splits the single interface—the 
juncture of built-in and free zone--into two interfaces; the juncture 10 of built- 
in and partially restrained zone, and the juncture 21 of partially restrained and 
free zone, as illustrated in Fig. 1. The question arises: By what factor is the 
bending stress reduced through the introduction of this buffer zone? We shall 
show that when the sand layer is regarded as having a spring constant K uni- 
form over its depth, a proper choice of x will result in a reduction of the maxi- 
mum bending stress to 2/9 of the value encountered when no transition zone is 
provided. If we also take the variation of x, from 0 at the surface to Kya, at 
the bottom, into account, then, naturally, the transition becomes much smooth- 
er, and by making Xm x very large (and having an adequate distance D avail- 
able) the bending stresses may be reduced to arbitrarily small values. Need- 
less to say, a similar buffer zone action can be provided not only by sand re- 
straint but also by the use of closely spaced metal reinforcing rings as indi- 
cated in Fig. 2. 

The contents of this paper are as follows: First we establish the formula 
for sand stiffness in terms of its more familiar characteristic, compressi- 
bility; we also give the formula for metal band stiffness. Then we write the 
expressions for the membrane behavior of the shell in each of the Zones 0, 1, 
2, considering the shell to be cut along the interface between each zone. Third, 
we determine the edge moments and edge shears developed at the interfaces 
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which maintain continuity between the shell sections in the three zones. 
Fourth, we calculate the stresses and deformations produced by the edge re- 
actions, and superpose these effects on the membrane results. The final re- 
sults are the curves drawn in Fig. 5 for radial displacement §, and in Fig. 6 
for meridional stress ox versus position nf, (counted from interface) for 
values 0, 1, 10, and ~ of the spring constant ratio 


x/k 
Here 


is the familiar stiffness expression of a shell of radius R and wall thickness 
h, and «is the soil stiffness, psi/in.; 8; is the attenuation length, as given by 
Equation (24), for Zone i; n is the number of attenuation lengths along the 
shell from the interface in question. 

The stresses plotted are given in units of 
(2) 
where op is the hoop stress in an unrestrained shell. The deformations plot- 


ted are given in units of Rp Jah 
pr (3) 


where bp is the outward radial displacement in an unrestrained shell (sealed 
at the two ends). Similar curves also have been calculated (but are not shown) 
for 0.1, 2, and 100, as well as for the resultant hoop stress, Og. For the Pois- 
son ratio, UW, the value 0.3 was used. 

The composite picture, Fig. 7, illustrates the main conclusions of the __ 
paper. Using subscripts i, 0 for inner and outer fiber, and subscripts 10, 21 
to indicate action arising in the vicinity of interfaces 10 and 21, Fig. 7 plots 
the highest values of %9210 2nd the lowest values of 


VETSuS P, as derived from Fig. 6 for and a similar 
(not shown) figure for og. All other o, and og stresses (such as the lowest 


value of oxj9j, or the highest value 0997; are intermediate to the upper and 
lower 0, curves, and to the upper and lower og curves, respectively. Thus 
the heavy contour in Fig. 7 indicates the dependence of the stress range on 
the stiffness ratio p. It is seen that the inner fiber meridional stress at the 
interface 10 is the largest stress for stiffness ratios from P= 0 to P= 2; in 
this range 0xjo;j decreases from 2.04 a, to 1.04 0, and the bending compo- 


nent of oxj{$j decreases from 1.54 Op to 0.54 Op. In the range p = 2 to P = 16 


the hoop stress in the outer fiber near interface 21 is largest; it has an es- 
sentially constant value; 0g3{, = 1.04 o). In the range P= 16 top = ~ the inner 


fiber meridional stress near interface 21 is the largest: 0,37; rises from 
1.04 dp to 2.04 op. 

We may refer to the range p= 2 to P= 16, where og = 1.04 % is the maxi- 
mum stress, as the optimum range, and to the point / = 3.8 where 0,j0j = 
0x3]; = 9.83 o, as the optimum point. For this value of pthe maximum bend- 
ing stress experienced by the shell is 0.33 o,, i.e., it is only about 1/4.5 of 
the maximum bending stress that develops when P= 0 or 


4 


Stiffness of Soil® 


Stiffness of soil is usually specified in the form of an e - p curve as illus- 
trated in Fig. 3a, derived from a compression test, illustrated in Fig. 3b. 
The soil sample is confined in a ring, and it is compressed by means of a 
rigid slab through the application of pressure p. Let e be the original void 
ratio of the soil, i.e. let 


where V is the total volume of the sample, Vy the total volume of voids. 

Fig. 3a plots a typical curve showing the decrease in void with an increase in 
pressure. For the curve shown, the void change in the range from p = 0 to 
150 psig is Ae = 0.05. If, as shown in Fig. 3b, the occupied volume is per- 
mitted to change only through a change in a single dimension, L, then we may 
write 


r Ar r 
Consequently, the spring constant for the arrangement shown in Fig. 3b is 
= Ac (6) 


If the depth D of the sand layer in Fig. 1 is three or more times the width 
L, then the existence of a transition layer below the interface 21, about L 
deep, where the soil stiffness gradually increases from zero to its full value, 
x, may be ignored (this is conservative). If the width L is small also in com- 
parison with the shell radius R, then we may regard compression of the sand 
as a one-dimensional affair. It follows that for e = 0.8, Ap = 150 psi, Ae = 
0.05, and a width L = 2 ft. of the sand layer (and assuming also D >3L, R>> 
L), the soil provides a stiffness int 


1.8 = 


«© 226 fre! Sin (7a) 


Consider furthermore a pressure vessel of radius R = 20 ft. wall thick- 
ness h = 1 in, and modulus of elasticity E = 30 x 106 psi. Its stiffness is 


ka = « Sec (2 =! (7b) 
Thus, for the foregoing numerical values, 520 psi pressure will expand the 
vessel 1 in. and 225 psi pressure will compress the sand layer 1 in. Con- 
sequently, ratio of soil stiffness to pressure vessel stiffness is 


6520 2«© 0.43 (8) 


Of course, one must be careful not to take the stiffness Equation (6)—which 
ignores nonlinearity, surface behavior, uneven compacting, pockets, change 
with atmospheric conditions, water-logging, etc.—too literally. If we wish to 
have a ratio p = 4.3, we can not, in practice, insure it by reducing L from 24 
in. to 2.4 in., because, in pouring sand into such a narrow channel, we may 
no longer ignore the departures from assumed uniform behavior. But for a 
layer 2 ft. wide, and more than 4 ft. deep, we are justified in concluding from 
Fig. 7, abscissa p = 0.43, that the maximum stress in the container is 


5. This section is based on Article 13 in “Soil Mechanics,” by K. Terzaghi 
and R. B. Peck, Wiley and Sons, 1948. Fig. 3a represents conditions for 
dense sand. 
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(9) 


of which 0.50 Op is meridional tensile stress caused by the pressure, and 
1.1 Gp is bending stress caused by the sudden change in restraint at the inter - 
face 10 where the concrete and the sand layer meet. 


Stiffness of Metal Bands 


The stiffness provided by thin-walled metal rings is well-known, 
where h, is the ring wall-thickness, R, is the mean ring radius, and y is the 
ring height t divided by the ring spacing s, as indicated in Fig. 2. 


Membrane State 


We consider the shell dissected at the interfaces 01 and 12, and let each 
section acquire its own membrane deformations and stresses under the in- 
ward pointing normal load Z, psi, as indicated in Fig. 4a. For Zone 1 this is 


and leads to the hoop force per unit length of cylinder 
In the above ¢ is the radial outward displacement, 0, the hoop stress, and Ng 
the hoop force. On the other hand, the meridional force and stress depend 


only on p, 
= 
| (13) 


Here op introduced earlier in Equation (2), is the pressure hoop stress in an 
unrestrained shell. 
From the expression 
= 3, kK = \ N ~ aN (14) 


4 + / 
of the hoop strain we eliminate § by means of Equations (11), (12), and (13), 
and obtain 


NA A op Lie g) (15) 
From Equations (13) and (14) it now follows that 


R p/ Ea) 


where dp, introduced earlier in Equation (3), is the radial displacement in an 
unrestrained shell. Placing P = 0 in these expressions we obtain the familiar 
formulas 


@ 


which pertain to the unrestrained zone; placing P= ~ we obtain 


which pertain to the fully restrained zone. We assumed that the concrete 
foundation prohibits radial expansion but does not prevent longitudinal stretch- 


ing. 


2#t66 
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Edge Moments and Forces which Eliminate Misalignments 


Henceforth we shall use subscripts7 for the membrane results derived in 
the previous section, to distinguish them from the effects of edge load actions 
to be introduced in the following: In Fig. 4 we show the misalignments 


* 
594 = - Sim = +p) 


= Som = 6, +p) 


in displacement which arise at interfaces 21 and 10 due to membrane action 
only. Since the membrane state gives rise to no rotations 


X0m = Xim = Xam = 9 (20) 


(19a , b) 


there are not rotational misalignments, 


X91 = Xam + Xm = 9 


(21) 


* 
= Xim + = 


Note that we consider a rotation X positive if it rotates inward as viewed 
from the common edge, as shown in Fig. 4b. Likewise we regard shear force 
H as positive if it points outward as viewed from the common edge. On the 
other hand we define positive bending moment M as inward, and positive dis- 
placement 6 as outward, irrespective of the viewer's location. As a conse- 
quence the displacement misalignments in Equation (19) involve the difference 
of the membrane displacements, the rotational misalignments in Equation (21) 
involve the sum of the membrane rotations. The foregoing is in accordance 
with the sign conventions introduced by Horvay and Clausen.® 

It is readily seen that to eliminate the misalignments, edge shears H 2» 

iw 1 14, 9 and edge moments 2, 14,74 must be developed 
(1 u = upper edge of Zone 1, 1! = lower edge of Zone 1) which are in equilib- 


rium 
M2 H2+Hiy=09 


(22) 
M 
and which give rise to deformations at the upper edge of Zone 1 
= (-A1u +H X1u = iu - 5H 1u41) (23) 


Similar expressions 6,9, 4); 59, 2; 49, Xo also hold at the lower edge 

of Zone 1, at the lower edge of Zone 2, and at the upper edge of Zone 0, res- 

pectively. The bar above the letter denotes the “edge value” of deformations 

6 and x ;a change from Latin to script letter signifies the edge value of 

moment and force My and H. 
In the foregoing 


1/4 
= (4/K) / (24a) 


is the “attenuation length” of Zone 1 of the shell, where 


6. “Membrane and bending analysis of axisymmetrically loaded axisymmetri- 
cal shells,” by G. Horvay and I. M. Clausen. J.A.M., in press. 
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= 


K=k+k (24b) 
denotes the combined spring restraint provided by the shell wall and the soil, 
= En3/12 (1 - (24c) 

is the bending stiffness of the shell, and 


fy, (24d) 


is the “flexibility” of the shell. 
For shell Zones 2 and 0 the corresponding quantities are obtained by sub- 
stituting into Equation (24b) values of k = 0 and K= ~, respectively. 
From the four continuity requirements 
- $2 = 991, Xiu + =0 as) 
- = Xo + = 0 
and using the four equilibrium Equations (22), we solve for the 8 unknown edge 
reactions, and find in terms of the ratio 
/ /4 


B = (Ky /K,)" = (1+ (26) 


the values 
2(1 - 


“(1 +B +B? 49) 
4p bp 
u 
(1+8+62+83)(1+p) 


a 2 6 


4 6 


Deformations and Stresses Produced by Edge Loading 


The edge actions /7/ w Wf iy Manifest themselves, in accordance with 
Table 2° in the radial depiscdtbent 


V2 sin (nyy - 7/4) + H Cos ny y| (28a) 
lu 


6 
lu 2 
and in the meridional bending moment 


-n 

Mx iu =e V2 cos (nyy - 7/4) - Huh sinn,,) (28b) 
in Zone 1, at a distance nj yf; from interface 21. Similar relations also hold 
for the quantities derived from //9,..., 7 g- Knowing 6 and My we obtain 
the meridional bending stress from 


Aas 
4 | | 
677-6 
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Oy = + 6 My/h? (29a) 


and the hoop stress from 


ug = E6/R + 6UM,/h? (29b) 


The upper sign refers to inner and the lower sign to outer fiber stress. 
Resultant Stresses and Displacements 


Referring to the effects of edge-load actions 7 9,..., H o listed in the 
previous section by subscript e, we determine the resultant deformation and 
stresses in Zone 1 of the shell taking 


6; dim + Die 
= "xim * (30) 
% 1 = * ole 


Corresponding quantities in the other zones are similarly formed. Expres- 
sions for de; Sxe, %ge were listed in the preceding section, and expressions 
for dm, %m » %m Were listed in the section on “Membrane State.” 

We plot the resultant deformation 6 and the resultant stress ox, so obtained 
throughout the shell, in Figs. 5 and 6 versus distance from interface, for p = 
0, 1, 10,. Our assumption that Zone 1 is several (at least 3) attenuation 
lengths long insures that what happens at one interface has no effect on the op- 
posite interface of this zone. 

From the curve of 0, shown on Fig. 6, and from a similar curve (not shown) 
for og, we derive the composite Fig. 7, plotting the maxima and minima of the 
respective stresses. For instance, we see from Fig. 6 that the maximum value 
which the meridional stress reaches in Zone 1 for Pp = 1 is 0x max * 1.2640, in 
the inner fiber at interface 10, and ox max = 9.660, in the inner fiber at ony 
distance from interface 21. These two values are plotted at p= 1 on the oxJpj 
and the 0,97; curves. 
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